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ABSTRACT 

The reionization of the Universe has profound effects on the way galaxies form and on 
their observed properties at later times. Of particular importance is the relative timing 
of the reionization history of a region and its halo assembly history, which can affect 
the nature of the first stars formed in that region, the properties and radial distribution 
of its stellar halo, globular cluster population and its satellite galaxies. We distinguish 
two basic cases for the reionization of a halo - internal reionization, whereby the stars 
forming in situ reionize their host galaxy, and external reionization, whereby the pro- 
genitor of a galaxy is reionized by external radiation before its own stars are able to 
form in sufhcient numbers. We use a set of large-scale radiative transfer and structure 
formation simulations, based on cosmologies derived from both WMAP 1-year and 
WMAP 3-year data, to evaluate the mean reionization redshifts and the probability 
of internal/external reionization for Local Group-like systems, galaxies in the field and 
central cD galaxies in clusters. We find that these probabilities are strongly dependent 
on the underlying cosmology and the efficiency of photon production, but also on the 
halo mass. There is a rapid transition between predominantly external and predom- 
inantly internal reionization at a mass scale of ~ IQ^^Mq (corresponding roughly to 
L* galaxies), with haloes less massive than this being reionized preferentially from 
distant sources. We provide a fit for the reionization redshift as a function of halo 
mass, which could be helpful to parameterize reionization in semi-analytical models 
of galaxy formation on cosmological scales. We find no statistical correlation between 
the reionization history of field galaxies and their environment. 

Key words: cosmology:theory - radiative transfer - large-scale structure of Universe 
- galaxies:evolution 



1 INTRODUCTION 

The lack of Gunn-Peterson troughs in the spectra of dis- 
tant galaxies and QSO's show that the Universe has been 
almost completely ionized ever since redshift z ~ 6. This 
reionization was a consequence of the ionizing radiation 
emitted by the first galaxies, which started forming at 

2 > 20 through the gravitational growth of the primor- 
dial density fluct uations. The recent 3-year results from the 
WMAP sateUite (jSpergel et al.ll2007f ) indicated that reion- 
ization was well-advanced by z ~ 10 and was thus fairly 
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extended in time. Reionization had some profound impli- 
cations on star and galaxy formation. First, stars form- 
ing in a reionized medium of low metallicity are probably 
considerably less massive than s tars forming in a neutra l 
medium of low metallicity (e.g. Ijohnson fc Bromml |2006| ) 
and follow a different evolutionary path. Second, the for- 
mation of low-mass galaxies is suppressed in the presence 
of an ionizing background. The shallow gravitational wells 
of the small haloes cannot hold onto t heir gas if the halo 
virial tem perature i s below Tyir ^ 10^ K (Dekel & Silk"l986|, 
.Shapiro etal]|2004 llliev. Shapiro, fc Raga.2005 ). while for 
larger haloes (10*K <^ Tvir <^ 10^ K) the gas accretion 
rate decreases, and the cooling time of the gas increases. 
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thereby suppressing further star formation (lEfstathioul 
19921: iThoul fc Weinber3ll996l: iNwarro fc Steinmetj|l997l : 



Bullock et al.ll2000l : iDiikstra et al.ll2004h . Therefore, reion- 
ization is a promising candidate for reconciling the discrep- 
ancy between the number of observed low luminosity satel- 
lite galaxies around the Local Group and the CDM simu- 
lation pre dictions, often refe r red to as the 'missing satellite 
problem' (Kl vpin et alll 19991 : iMoore et al.llT999l : ITuIIv et^ 
l2002i ). Third, reionization may have truncated the formation 
of the blue, old globular clusters that are still observed today 
l|Santoj|2003l :lM oore et al ] |2006l ). Since massive galaxies will 
be reionized considerably before reionization is completed 
in the entire IGM, it is interesting to ask when and how a 
galaxy of a given mass, located in a given environment, will 
be reionized. 

Generically, there are two ways in which a galaxy can 
be reionized. The process could be driven by the radia- 
tion emanating from its own stars forming in situ ('inter- 
nal reionization') or the galaxy or its progenitors could be 
overrun by the global ionization fronts created by objects 
which collapsed earlier ('external reionization'). The second 
scenario is expected to occur e.g. for low mass satellites 
that have been close to large galaxies like the Milky Way 
since before reionization, while the first scenario should be 
typical for very massive galaxies in the central regions of 
clusters. It is less clear what scenario should be expected 
for intermediate-mass galaxies. For example, a field galaxy 
which resides close to a massive cluster could have been 
reionized by the radiation coming from the proto-cluster be- 
fore the galaxy itself formed. The concept of a non-uniform, 
de nsity-dependent re i oniza tion redshift has been introduced 
by iBarkana fc Loebl ifioCMlV They point out that external 
reionization may become important in underdense regions 
and voids, thus causing the relatively low number of dwarf 
galaxies in these regions. 

In this work we study the reionization history of haloes 
and its dependence on the halo mass and environment. To 
this end, we explore data from series of large-scale struc- 
ture formation and radiative transfer simulations of cosmic 
reionization. We re-simulate the same initial conditions at 
lower resolution which allows for tagging of the particles and 
thus for following the detailed merger and accretion history 
of haloes from the reionization epoch to the present. This 
allows us to identify the progenitors of massive cluster galax- 
ies, field galaxies and Local Group-like systems and simulta- 
neously follow both their reionization history and their mass 
growth history. 

In this paper, we will use two difi'erent flat 
ACDM cosmologies, the first with parameters 
a}^,nA,ilb,hy j8, n) = (0.27,0.73,0.044,0.7,0.9,1) 
jSpergel et al.l |2003| . hereafter WMAPl) and 

rn^,nA,ni,,/i ff8,n ) = (0.24,0.76,0.042,0.73,0.74,0.95) 

jSpergel et ahlboOTl . hereafter WMAP3). Here t^m, Oa, and 
S7f, are the total matter, vacuum, and baryonic densities 
in units of the critical density, h is the Hubble constant 
in units of 100 km s~^Mpc~^, erg is the standard deviation 
of linear density fluctuations at present on the scale of 
8h~^Mpc, and n is the index of the primordial power 
spectrum. 

This paper is organized as follows: In § [5] we summa- 
rize our simulations. In § [3] we present our techniques used 
to select fleld haloes and follow the reionization histories of 



Table 1 . Overview of the three different codes and the resolutions 
used in our WMAP3 simulation (WMAPl is similar). Note that 
the spatial resolution in PMFAST is better by at least an order 
of magnitude as far as the position of halo centers are concerned. 
In the PKDGRAV simulation, the spatial resolution is not clearly 
defined and was assumed to be equal to the force resolution. 





code 


mass res. 


spatial res. 


high res 


PMFAST 


2 xlO^ /i-^Mq 


~ 30 h'^kpc 


low res 


PKDGRAV 


9.96 xlO** h-iMg 


~ 10 h-^kpc 


RT 


C^-ray 




~ 0.5 h-^Mpc 



individual haloes. In §|4]we present our results for the reion- 
ization histories of field galaxies, cD galaxies in the centers 
of galaxy clusters and Local Group-like objects. In § [5] and 
§ [6] we discuss potential caveats and discuss the implica- 
tions of external vs. internal reionization on the observable 
galaxy characteristics. Finally, in §(7] we give a summary of 
our results. 



2 SIMULATIONS 

Our simulations follow the evolution of a comoving simula- 
tion volume of (100 ^~^Mpc)''. Our basic metho dology and 



simulat i on parameters are desc r ibed in detail in 



Iliev et ahl 



(|2006bl ): [Meilema et al.1 (l2006bl ): llliev et al.1 l|2007h . Here we 
provide just a brief summary, in addition to a detailed dis- 
cussion of the elements which are first presented in this work. 
See Table[l]for listing of the codes used in this work and their 
effective mass and spatial resolutions. 



2.1 Structure Formation Simulations 

We start by performing a very large dark matter only sim- 
ulation of early structure formation, with 1624^ ~ 4.3 bil- 
lion particles and 3248^ grid cellfl using the code PMFAST 
erz et al1l2005l ). This allows reliable identification (with 
100 particles or more per halo) of all haloes with masses 
~ 2 X l(f ^"'^Mq or larger. We compile and save the halo 
catalogues, in up to 100 time slices starting from high red- 
shift (2 ~ 30) to the observed end of reionization at z ~ 6. 
In what follows, we will refer to this simulation as HR (high- 
resolution) simulation. This HR-simulation and in particular 
the halo catalogues based on it are used as the framework 
for our radiative transfer (RT) simulation which will be ex- 
plained in detail in section 12.21 In order to maximize the 
mass resolution, the dark matter particles are not labelled 
in the HR simulation, which does not allow us to identify 
the progenitors of any given structure at 2: = 0. We resolved 
this problem by re-simulating the same initial conditions at 
lower resolution (called LR simulation hereafter) which al- 
lows for tagging of the particles. This enables us to identify 



1 3248 = N^odes X (512 - 2 x 24), where N^^ds, = 7 (with 4 
processors each), 512 cells is the Fourier transform size and 24 
cells is the buffer zone needed for correct force matching on each 
side of the cube. 
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Figure 2. Reionization colour plot, showing an example of external reionization. Symbols/Colours have the same meaning as in Figure[T] 
At 2=7.346, the number of haloes is so large that we do not show them anymore in the plot. 



the progenitors of haloes at 2=0 and to simultaneously fol- 
low their reionization and mass growth histories. To obtain 
the LR simulation, the original HR PMFAST initial con- 
ditions (with 1624^ particles) were generated again using 
the same random seeds. They were then coarsened to 406^ 
particles by averaging the properties (positions and velocity 
components) of every 4^ neighboring particles. We evolved 
these coarsened initial conditions from z = 100 to z = us- 
ing PKDGRAV, a fast p arallel tree co de written by Joachim 
Stadel and Tom Quinn (jStadelll200ih . 



2.2 The Radiative Transfer Simulations 

We proceed by applying our radiative transfer code to the 
HR-simulation described above. Since radiative transfer sim- 
ulations at the full grid size of our N-body simulations are 
still impractical, we follow the radiative transfer on coarser 
grids, of sizes 203^. All identified haloes are assumed to be 
sources of ionizing radiation and each is assigned a photon 
emissivity proportional to its total mass, M, according to 



iV^ = A 



A-mt 



(1) 



fimpt sQo 

where ts is the source lifetime, is the proton mass, fi is 
the mean molecular weight and f.y is a photon production ef- 
ficiency which includes the number of photons produced per 
stellar atom, the star formation efficiency (i.e. what frac- 
tion of the baryons are converted into stars) and the escape 



fraction (i.e. how many of the produced ionizing photons es- 
cape the haloes and are available to ionize the IGM). Since 
we do include in the halo catalogue only haloes defined by 
more than 100 particles, this criterion introduces a mini- 
mum mass for haloes in order to act as ionizing sources. For 
our choice of the simulation parameters this corresponds to 
a mass threshold of Mth = 2.5 x lO" h'^MQ (WMAPl) 
and Mth = 2.2 x 10^ /i~^M0 (WMAP3). Note that this 
mass corresponds roughly to the halo mass above which a 
galaxy forming i n it is largely unaffected by reionization (see 
llliev et aUboOTl . for discussion and references). 

The radiative transfer is followed using our fast and 
accurate ray-tracing photoioniza tion and n on-equilibrium 
chemistry code C'^-Ray ( McUcm a et al.l bOO Gal . The code 
has been tested in detail for correctness and accuracy against 
available analytical solutions and a number of other cos- 
mological radiat ive transfer codes (|Mellema et al ] l2006al : 
llliev et al.ir2006ah . The radiation is traced from every source 
on the grid to every cell using short-characteristic ray- 
tracing. The runs and n o tation a re the same as t he one s 
presented ir i llliev et al.l l|2006bD : iMellema et al] (|2006bD : 
llliev et al] (|2007 ). These include simulations with both 
WMAPl and WMAP3 background cosmologies. The runs 
for a given cosmology share the same underlying N- 
body simulation (but with different random realizations for 
WMAPl and WMAP3), and each adopts different assump- 
tions about the source efficiencies and the sub-grid level of 
gas density fluctuations. In part of the simulations, a sub- 
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Figure 1. Reionization colour plot, showing an example of in- 
ternal reionization. Redshifts are denoted on the top of the plots. 
The reionization status of the gas is shown in colours, dark blue 
showing completely neutral gas, orange completely reionized gas. 
Black rectangles denote the position of identified haloes, the white 
dashed circle denotes the position of the LG progenitor, identi- 
fied by its densest particles. The ionization state of the gas was 
identified in a slice with a thickness of one cell of the simulation 
outputs at different redshifts, with the position of the slice cen- 
tered around the Local Group progenitor in consideration. Iden- 
tified haloes are shown if they are within a box of volume 20^ 
cells centered around the Local Group progenitor (i.e. can be up 
to 10 cells offset from the slice shown in the pictures). 



grid gas clumping is included, C{z) 
evolves with redshift according to 

/ \ r.'v /ICC -0.1142+0.001328 

Csubgridl^j = 27.466e 

in WMAPl cosmology and as 



Csub-grid(^) = 26.2917e' 



-0.1822Z+0. 003505 z^ 



(n2)/(n)2, which 
(2) 

(3) 



for WMAP3 cosmology. These fits were obtained from an- 
other set of two high-resolution PMFAST N-body simula- 
tion, with box sizes (3.5 h^^ Mpc)^ and a computational 
mesh and number of particles of 3248^ and 1624"^, respec- 
tively. The inclusion of gas clumping results in a slower prop- 
agation of the reionization fronts and a delay of the final 
overlap. 

In what follows, we will concentrate on three cases, 
two 'extreme' models and an intermediate one. The extreme 
cases (in terms of Zov, the redshift at which the final over- 
lap is reached), are the following: WMAPl-f2000 (WMAPl, 
= 2000, no sub-grid gas clumping) and WMAP3-f250C 
(WMAP3, fj = 250, with sub-grid gas clumping). These 
two simulations should roughly bracket the allowed range 
in cosmological and reionization parameters. Among all of 
our radiative transfer simulations, the WMAPl-f2000 sim- 
ulation has the earliest redshift of overlap, Zov^H-S, while 
simulation WMAP3-f250C has the latest overlap, Zov ~ 6.5. 
We add a simulation in between those two cases: WMAPl- 
f250C (WMAPl, f-y = 250, with sub-grid gas clumping), 
resulting in 2:ov=8.2. This allows us to study the effect of 
changing only the cosmology or changing only the reion- 
izing photon production efficiency. The lower efficiency of 
/t = 250 appears more realistic at present, and the corre- 
sponding cases match the late end of reionization at 2 ~ 7 
well. We note that the delayed structure formation due to 
the lower power spectrum normalization, tilt and lower mat- 
ter density in the favoured WMAP3 cosmology naturally 
predicts the lower integrated Thomson scattering optical 
depth the same measurements found, and thus do not seem 
to require a lower r eionizing phot on production efficiency 
ll Alvarez et all l2006l. but see also iHaiman fc Brvanl I2OO6I : 
IShuU fc Venkatesanll2007l ). 

For the WMAPl simulations, transmissive boundary 
conditions on the ray-tracing were used. As a consequence, 
haloes close to the boundaries of our computational volume 
will be somewhat less affected by external radiation, partic- 
ularly towards the end of reionization. Thus, in order to be 
conservative, we do not include in our analysis any haloes 
that are within 10 /i~^Mpc from the simulation boundary. 
The WMAP3 simulation utilized periodic boundary condi- 
tions, which allows for all haloes to be used in our analysis. 



3 METHODS 

3.1 Matching the HR and LR Simulations 

In order to investigate the reionization and formation his- 
tory of a given halo, we first identify it at z—Q in our LR 
simulation, using a Spherical Overdensity Halo Finder (see 
iMaccio et al.|[2007l for a detailed description). We then track 
the particles of the halo in consideration back to higher red- 
shift in the LR simulation and match them to the HR sim- 
ulation in order to analyze its formation and reionization 
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history. This can only be done if the density distribution in 
the HR and LR simulation agree well; since two quite diflFer- 
ent methods and resolutions are used for them, this merits 
checking and evaluating the closeness of the match between 
the two. If they match well, the location of the haloes which 
can be identified in the LR simulation should be similar to 
the location of the most massive haloes in the HR simulation 
(where the total number of identified haloes at any redshift 
is much higher than in the LR simulation). As a test, we 
identify for each halo in the LR simulation which is found 
by our Spherical Overdensity Halo Finder the nearest halo 
in the HR simulation which has a similar mass (within a fac- 
tor of 2). At z=6 in the WMAP3 simulation, we find that 
the median distance between the haloes in the two cata- 
logues is 0.3 Mpc , with a similar result for the WMAPl 
simulation. Since this distance is below the grid size of the 
RT code, we consider this reasonable agreement for our cur- 
rent purposes. However, for around 3 percent of all cases, 
no nearby halo with similar mass is identified. This is not 
surprising given the fact that we use two difi'erent halo find- 
ers and have a large difference in resolution. Note that the 
median distance and the number of unmatched haloes de- 
creases with increasing redshift. This reflects the fact that 
the DM distribution in the HR and LR simulation becomes 
more similar at high z, since the effects of different force 
resolution, etc., are diminished. As all of our comparisons 
between LR and HR simulations are done at times when 
the Universe is not yet completely reionized, i.e. 2 > 6, we 
consider the median distance and the fraction of unmatched 
haloes quoted above an upper limit. 

In what follows, we will explain in more detail how the 
formation and the reionization history of a given z — halo 
is investigated. We mark the particles that will end up in 
a given halo at z = and track them back to simulation 
outputs at higher redshifts. At any redshift, we identify the 
center of the halo progenitor as the densest particle (smooth- 
ing over the 32 closest neighbours) in the marked set and 
measure the fraction of reionized gas in the central cell and 
the 26 neighbouring cells, using the radiative transfer grid 
consisting of 203^ cells. In this way, an offset of one radiative 
transfer (RT) cell (corresponding to ~ 0.5/i~^Mpc) in each 
dimension is allowed between the center of reionization in 
the HR-simulation, and the location of the densest particle 
in the LR-simulation. If the ionized fraction in any of those 
cells is higher than a threshold value of Xth = 0.7, we de- 
fine the halo as reionized and the corresponding reionization 
redshift Zr is set equal to the current output redshift. The 
particular value of xth has a negligible impact on the final 
results thanks to the sharp transition across the ionization 
fronts (see Fig. [T]and[2|. 

Next we check at which redshift the most massive pro- 
genitor of the halo in consideration first reaches a halo mass 
of Mth (as defined in § 12. ip . and define this as redshift of 
formation, Zf. To do this, we again trace back the halo pro- 
genitor particles and identify the densest particle at each 
redshift. Again allowing for an offset of one cell, we consider 
a halo as formed if it is identified by our halo finder based 
on the HR simulation. Only for less than 2 percent of all 
haloes, no formation redshift could be identified. They were 
removed from the sample. 

Ionizing radiation in the ray tracing simulation can, by 
definition, only be produced by haloes that have been iden- 



tified; this allows us to clearly separate internally reionized 
haloes from externally reionized ones. If Zr > Zf {zr < Zf), 
the halo has been reionized externally (internally) . For illus- 
trative purposes, in Figures[T]and[2]we show the evolution of 
the ionization fronts around an internally and an externally 
reionized Local Group-like object, respectively. 

Our method can occasionally fail if the halo progenitor 
under consideration has a very large extent (i.e. is larger 
than two cells) and has two dense centers. If one of those 
is identified as densest particle by our smoothing algorithm, 
while the other is recognized by the high-resolution halo 
finder, external reionization will be always assumed. Eye in- 
spection indicates that this happens in less than roughly 
10% percent of all cases we have identified as external reion- 
ization. On the other hand, if the reionization and the for- 
mation of a halo occur in-between two timesteps, we always 
assume internal reionization. Although external reionization 
is also possible it is unlikely. We did not find a way to solve 
these two issues, but we expect that their impact on our 
results should be modest, since these effects will both only 
occur for a small fraction of our haloes. 

3.2 Sample Selection 

We have identified dark matter haloes at z—0 in our LR 
simulation using our Spherical Overdensity Halo Finder. 
From this halo catalogue we selected three sets of objects 
for further consideration: 1) field haloes, 2) central objects 
of massive clusters, likely to correspond to central cD galax- 
ies and 3) Local Group-like systems. As a subset of the field 
haloes, we also selected a sample of haloes likely to host an 
L* galaxy. 

The sample sets of such systems were chosen as follows. 
The field halo samples consist of all the haloes with masses 
between 1.8 x lO" H'^Mq and 3.7 x lO" /i^^M© that have 
been found by the Spherical Overdensity Halo Finder in both 
simulations. Such haloes are always isolated in the sense that 
they are not subhaloes in a larger group or cluster. In the 
WMAPl-f2000 and the WMAPl-f250C simulation, haloes 
with distances closer than 10 Mpc from the edges of the 
simulation box were removed from the sample. In this way, 
we have identified 9683 haloes in the WMAP3-f250C simula- 
tion, and 7346 haloes in the WMAPl-f2000/WMAPl-f250C 
simulations. A sample of haloes lik ely to host an L* galaxy 
was chosen as follows. According to lMadgwick et aU ( 2002t ). 
L* - galaxies have a 6j-luminosity of 1.1 x 10^° h~^LQ . A cen- 
tral or isolated L* galaxy is therefore most likely to be found 
in a halo with log(M) = 12.2 [ /I'^Mo] in a WMAP3, and in 
a halo with mass log(M)=12.3 [ H'^Mq] in a WMAPl cos- 
m ology, following the re l ation between Lcen and halo mass 
in Ivan den Bosch et al.1 (|2007| ). The difference in the two 
cosmologies is mainly due to the different halo clustering 
properties. We select for our L* galaxy sample haloes with 
mass log(M)=12.1-12.3 [ /i~^Mq] for the simulation based 
on WMAP3, and with log(M) = 12. 2-12.4 [ fe'^M©] for the 
simulations based on WMAPl. We find 437 such haloes in 
the two simulations based on the WMAPl cosmology, and 
575 haloes in the WMAP3-f250C simulation. 

The second set of haloes consists of the central objects 
in the 20 most massive galaxy clusters in each simulation. 
These massive objects are likely to correspond to central 
cluster cD galaxies. Finally, the third halo sample consists 
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Table 2. Summary of the results. Listed are the mean redshift 
of formation, Zf, the mean starting redshift of the halo reioniza- 
tion, Zr, and the mean redshift at which 70 percent of the final 
halo mass is reionized, 270%; the fraction of objects reionized ex- 
ternally, /ext, the average time interval between formation and 
reionization, At = tf — tr (in Myr), and At*, which is the same 
time difference in units of the Hubble time when 50 % of the mass 
in the universe is ionized, for all of our simulations, as labelled. If 
At is negative, haloes are on average reionized before formation. 





2/ 


Zr 


^70% 


/cxt 


At 


At* 


WMAPl-f2000 


field haloes 


10.8 


13.4 


13.1 


0.79 


-114 


-0.237 


L* halo sample 


12.6 


13.6 


12.85 


0.62 


-36 


-0.075 


central cDs 


16.7 


15.5 


14.1 


0.0 


26 


0.055 


LG sample 


12.8 


13.3 


12.9 


0.44 


-21 


-0.044 


WMAPl-f250C 


field haloes 


11.0 


11.7 


10.4 


0.4 


-33 


-0.053 


L* halo sample 


12.6 


12.4 


10.9 


0.17 


8 


0.013 


central cDs 


16.7 


15.5 


11.8 


0.0 


26 


0.042 


LG sample 


12.8 


12.4 


10.8 


0.1 


13 


0.021 


WMAP3-f250C 


field haloes 


8.2 


8.9 


8.2 


0.56 


-60 


-0.062 


L* halo sample 


9.3 


9.4 


8.6 


0.27 


-9 


-0.009 


central cDs 


13.0 


12.6 


9.6 


0.0 


15 


0.016 


LG sample 


9.1 


9.2 


8.3 


0.24 


-8 


-0.008 



of all Local Group-like systems in our volume. These are 
identified as binary systems of dark matter haloes satisfy- 
ing th e following criteria (see lMaccio. Governato fc Horelloul 
I2OO5I . for details): 

(i) The mass of each of the two haloes is between 4 x 10^^ 
and 5 x 10^^ /i'^Mq. 

(ii) The distance between the two haloes is below 700 
kpc. 

(iii) The two haloes are moving towards each other. 

(iv) There is no massive neighbour with a mass above 
5 X 10" closer than 2.1 Mpc. 



4 RESULTS 

4.1 Reionization Histories of Galaxies of Different 
Mass and Environment 

On average, reionization predates formation for typical field 
haloes in all three simulation cases that we have considered. 
We compare the redshift of formation, Zf, the starting red- 
shift of the halo reionization, Zr, and the redshift at which 
70 percent of the final halo mass is reionized, Zro%, for all 
our halo samples (field haloes, haloes likely to host an L* 
galaxy, haloes likely to host a cD galaxy, LG like haloes) in 
Table (2] We note that by definition Zr > Zm% , and espe- 
cially for massive haloes these redshifts could differ signifi- 
cantly since the material that will end up in the final halo is 
generally spread out over a very large region in space, and 
therefore a fairly long time could pass from start to finish 
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Figure 3. The average formation and reionization redshift (calcu- 
lated converting redshift to lookback time and back) as a function 
of the halo mass for the WMAP 1-12000 (top panel), the WMAPl- 
f250C (middle panel) and the WMAP3-f250C (bottom panel) sim- 
ulation. Error bars/dotted lines denote the 68 percent confidence 
level. The black solid line is a fit to the average redshift of reion- 
ization, as explained in the text. Note that formation redshift here 
means the redshift at which M > M^^. 



of the reionization of that halo's material. As mentioned be- 
fore, Zf is the redshift at which the largest progenitor of the 
halo reaches Mth- 

We find that external reionization occurs for roughly 
80% of all field haloes in the WMAPl-f2000 simulation, for 
40 % in the WMAPl-f250C simulation and for ~ 55% in 
the WMAP3-f250C simulation. Apparently, our "interme- 
diate case" with a redshift of overlap in between the two 
other cases, has the lowest number of externally reionized 
haloes. This is caused by an interplay of two different ef- 
fects. First, the higher the reionizing photon production, 
the more effective is a halo in reionizing its surroundings. 
If /-,=2000, each halo quickly reionizes a large volume, con- 
taining a mass ~ 2000 times the halo mass, as recombina- 
tions are not very important. This explains the very high 
fraction of external reionization in the WMAPl-f2000 case, 
which indicates that in this case, few early collapsing haloes 
are responsible for most of the reionization of the matter 
in the universe. The influence of the cosmo logical frame- 
work is more complicated. As pointed out by I Alvarez et al.l 
(2006), structure formation is delayed by a factor of ~ 1.4 in 
(1+z) in a WMAP3 cosmology compared to a WMAPl cos- 
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mology. This means that structure formation becomes more 
extended in time in a WMAP3 cosmology, i.e. the time of 
formation of low and high mass haloes are further apart. In 
order for external reionization to occur, the radiation of an 
early collapsing halo has to reach other haloes before they 
have started to form. This window in time will be shorter 
in a WMAPl cosmology than a WMAP3 cosmology, which 
explains the lower external fraction in the former case. This 
means that for a WMAP3 cosmology and a higher reionizing 
photon production efficiency, the external fraction would be 
even above 80 %. 

Another important outcome of our analysis is that dif- 
ferences in formation time cannot be directly translated into 
a difference in reionization time (cfr. Table [2]). Due to the 
contribution of external reionization, field galaxies "catch 
up" with central cD galaxies when it comes to reioniza- 
tion. The time elapsed between formation and reionization 
(At) depends both on the cosmological model and on the 
reionization parameters. In the early-reionization scenario 
WMAPl-f2000, At for field haloes and LG-like objects is al- 
ways negative and maximal in absolute value, which means 
that these objects form a fairly long time after they have 
been reionized. Field haloes are reionized on average notice- 
ably earlier in WMAPl-f2000 model than in the extended 
reionization scenarios. This is to be expected considering the 
much faster and earlier reionization in this case, allowing 
for significantly fewer field haloes to form before reioniza- 
tion. Some L* haloes in the WMAPl-f2000 simulation also 
do not manage to reach a mass of Mth before their reion- 
ization. Thus the coeval formation and reionization of L* 
haloes observed in the other two simulation cases is broken 
for early reionization. WMAPl-f250C and WMAP3-f250C 
share the same reionization parameterization. Thus, differ- 
ences in At are due to the different cosmological models. The 
characteristic time scale for structure formation is roughly 
the current Hubble time, tniz). So, in an attempt to scale- 
out the cosmological dependence of At, we divided it by 
the Hubble time tu,5o% when the ionized fraction (in mass) 
is 50% to obtain At*=At/tH.5o%. Comparing the extended 
reionization scenarios, WMAPl-f250C and WMAP3-f250C, 
we see that At* is very similar for all objects. 

Haloes likely to host massive cD galaxies in centers 
of clusters are always internally reionized in all the differ- 
ent simulations. This reflects the strong inside-out nature 
of reionization. Reionization always starts from the highest 
density peaks, inside which cD haloes reside, and thus these 
are reionized early, regardless of the details of the reioniza- 
tion parameters and the background cosmology. 

In summary, the timing of reionization compared to for- 
mation for L* and field haloes is sensitive to the reioniza- 
tion parameters (ionizing efficiency of the sources and gas 
clumping) and rather insensitive to the background cosmol- 
ogy, while in contrast the reionization history of central cD 
galaxies is largely independent of either the reionization pa- 
rameters, or the cosmology. 

In Fig. 121 we show the mean values of Zr and 2/ as a 
function of the halo mass at z = 0. Averages are calculated 
by converting redshift to lookback time, averaging and con- 
verting the result again back to redshift. The reionization 
redshift as a function of halo mass at z=Q can be fitted by 
a simple polynomial: 



Table 3. Fitting parameters for (zr) (M), with {zr) = ax + bx^ + 
cx^ , where x=\o^{M / Mq) 





a 


b 


c 


WMAPl-f2000 


5.079 


-0.599 


0.022 


WMAP3-f250C 


0.187 


0.079 


-0.001 


WMAP3-f250C 


1.592 


-0.167 


0.008 




11.5 12 12.5 13 13.5 
log.oM [h-1 MJ 

Figure 4. The fraction of externally ionized haloes as a function 
of halo mass respectively for the three simulations considered. 

{zr) — ax + bx^ + cx'^ (4) 

where a;=log(M/M0), expressed without . We list the 
relevant values for a, 6, c in Table |3] 

There are clear trends for the redshift of formation 
and reionization as a function of halo mass. For lower mass 
haloes, reionization typically happens prior to their forma- 
tion (external reionization), but for higher mass haloes, for- 
mation tends to predat e reionization (inte rnal reionization). 
As mentioned above, lAlvarez et al.l ()2006h found that struc- 
ture formation is delayed by a factor of 1.4 in (l+z) in a 
WMAP3 cosmology. We roughly find a similar change in 
the average formation time of a halo, comparing WMAPl 
and WMAP3 cosmologies. The 68 percent confidence level 
around average formation redshift is considerably larger for 
a WMAPl cosmology than for a WMAP3 cosmology at low 
masses. However, this simply reflects the fact that a simi- 
lar interval in time corresponds to a much larger interval in 
redshift at higher redshifts. Measured in lookback time, the 
scatter around the average formation redshift is very similar 
in the two cosmologies and even slightly higher for WMAP3 
at low masses. 

In Fig. |4j we show the fraction of externally reionized 
haloes as a function of halo mass. It is apparent that the 
probability of external reionization is strongly dependent on 
mass, mainly because the formation redshift and the mass 
of a halo are correlated, and haloes which form later have 
a higher probability of being externally reionized. We note 
that even if on average a halo forms after it vicinity has been 
reionized, the fraction of externally reionized haloes can still 
be rather low. The reason for this is that for many internally 
reionized haloes, the redshift of formation and the redshift 
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Figure 5. Top Panel: The fraction of externally ionized haloes 
as a function of the overdensity by measuring only the mass in 
collapsed haloes with masses above 5 X lO^'^. The overdensity is 
calculated with respect to the mean matter density in the simu- 
lation box. Results for all three simulations are shown. Errorbars 
denote the l-a error. Bottom Panel: Same as shown top panel, for 
the WMAP3-f250C simulation only, but for three different bins 
regarding the mass of the halo in consideration. The way the sur- 
rounding overdensity is measured remains unchanged. For clarity, 
no errorbars are shown. 
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Figure 6. The average formation and reionization redshift (cal- 
culated converting redshift to lookback time and back) as a func- 
tion of the overdensity measured only for the mass in collapsed 
haloes with masses above 5 xlO^^, for the WMAPl-f2000 (top 
panel), WMAPl-f250C and the WMAP3-f250C (bottom panel) 
simulation. Dotted lines/errorbars denote the l-a error. 



of reionization are the same or at least very close. Thus, even 
a relatively small fraction of externally reionized haloes can 
shift the mean reionization redshift so that it predates the 
formation redshift. 

Generally, one might expect that the probability of ex- 
ternal reionization of a field galaxy could also depend upon 
the local overdensity around it, since the latter is closely re- 
lated to the abundance of sources of ionizing radiation near 
the field galaxy in consideration. In Fig. [5] (upper panel) 
we plot the fraction of externally reionized field haloes as a 
function of the surrounding overdensity in collapsed haloes 
with mass M > 5 x 10^^ /i~^Mq at z = 0. The overdensity 
is measured within a sphere of radius 10 h~^Mpc centered 
on the halo in consideration, does not include this central 
halo, and is measured relative to the average density of the 
simulation box. Contrary to what one might expect, there 
is no clear trend. 

This result persists if we change the mass limit on the 
surrounding collapsed haloes, or the radius of the sphere. 
In Fig. [51 bottom panel, we show for the example of the 
WMAP3-f250C simulation that the constant behaviour of 



the externally reionized halo fraction also persists if we di- 
vide the halo sample into different mass bins. As evident in 
Fig. [S] the reason for the constant behaviour of the external 
fraction with halo mass is that not only the reionization, but 
also the formation of a halo occurs at slightly lower redshifts 
in under-dense regions, since the local "clock" of structure 
formation is slower in voids. Therefore, we conclude that the 
environment of haloes does not change the way in which they 
are reionized, at least if the dark matter halo is isolated. For 
cluster and group galaxies being part of a larger dark matter 
halo, this statement might not be correct anymore. Satellite 
galaxies in massive clusters that are close to the center cD 
galaxy will most probably be reionized by the time when 
70 percent of the dark matter mass surrounding the cen- 
tral cD galaxy itself is reionized. This is clearly earlier than 
the typical redshift where galaxies in the field are reionized. 
However, the difference is probably less dramatic than the 
difference in formation redshift. One could speculate that 
this means that galaxies in clusters are reionized internally 
more often than galaxies in the field. We have found this 
to be the case at least for the resolvable subhaloes in clus- 
ters, which we have identified with our subhalo finder SKID 



The Local Reionization History 9 



_l 1 1 1 1 1 1 1 1 1 1 1 

- WMAPl-f250C 
10 % external 

: r h 


1 1 1 1 1 1 1 1 1 

external - 
seir 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 
" J " 

, , 1 Uv 1 . . 1 , . . 1 . " 


- WMAP3-ta50C 
24 % external 

■, , 1 , , , 1 1 — 




, 1 , , , 1 , — 1 ' 


, . . . , 


1 1 1 1 . 
1 ,■ 



16 14 12 10 8 6 16 14 12 10 8 6 



Figure 7. Comparison between formation redshift (2:^-) and reion- 
ization rcdshift (zr) for LGs that are externally (diagonal texture) 
and internally (grey area) reionized. Results are shown for the 
WMAPl-f250C and the WMAP3-f250C simulation. Again, 2^ is 
defined as the redshift where the mass of the halo exceeds Aith- 



(|Stade]|l200ll ). However, since they are not a fair tracer of the 
subhalo population, we do not discuss our results in detail. 



4.2 The Reionization History of the Local Group 

In this section, we study the reionization history of Local 
Group like objects. As binary systems. Local Group like ob- 
jects might be expected to have a different average forma- 
tion and reionization history than typical field galaxies. Our 
Local Group of galaxies also provides the most detailed set 
of observational data for testing the consequences of reion- 
ization on subsequent galaxy formation. We apply the same 
definitions of internal and external reionization we described 
above to our LG candidate sample. For WMAP3-f250C we 
find that 24% (32 out of 133 LGs) are externally reionized. 
For WMAPl-f250C, the probability of external reionization 
is only ~ 10% (6 out of 62 LG). As for L* and field galax- 
ies, WMAPl-f2000 shows the highest fraction of external 
reionized LGs: 44% (27 out of 62). In Fig. we show the 
distribution of the reionization and formation redshifts of 
LG candidates for the two simulations with the lower reion- 
ization efiiciency: WMAPl-f250C and WMAP3-f250C. As 
in Fig.|3]at the lower mass end, it is evident that the scatter 
in formation redshift is larger than the scatter in reioniza- 
tion redshift. This is caused by the contribution of external 
reionization. 

To summarize, the probability for LG-like objects to be 
externally reionized is lower than the corresponding proba- 
bility for low-mass field galaxies, which can be as high as 60- 
80%. Compared to typical L* galaxies, LG-like systems form 
slightly earlier, get reionized slightly earlier and are some- 
what less likely to be externally reionized than L* galaxies. 



5 POTENTIAL CAVEATS 

All of our simulations resolve only haloes with masses above 
Mth ~ 2 X 10^ M0 and these are the only sources which con- 
tribute to the reionization of the universe in this set of sim- 
ulations. However, in hierarchical structure formation sce- 
narios the smallest haloes form first and merge up. Thus it 
is expected that the first sites of star formation were "mini- 
haloes" with a mass roughly between 10^ and 10* h~^MQ. 
By definition the minihaloes have virial temperatures be- 
low 10* K, where the radiative gas cooling through hydro- 
gen and helium atomic lines is not efficient. Thus, these 
haloes have to rely on the much less efficient molecular 
line cooling for their star formation. However, the hydro- 
gen molecules are easily destroy ed by radiative feedback 
(e.g. iHaimanT Abel, fc Reell2000l ) and thus it is generally 
expected that most minihalo es were 'sterilized' and in most 
cases did not form stars (e.g. IShapiro et al. l l2oo4 lii iev et al.l 
[2005 . and references therein). Haloes with a mass between 
10* and 10^ Mq, on the other hand, which are also not re- 
solved in our simulations, are able to cool efficiently and 
may actua l ly be an important source of ionizing radiation. 
Illiev et al] (|2007l ) have performed a number of smaller-scale 
high-resolution simulations of reionization, with a simulation 
volume (35h~^Mpc)^, which include radiation from these 
smaller haloes. The most important effect is that reioniza- 
tion starts much earlier than without those sources. How- 
ever, due to Jeans-mass filtering in the ionized regions and 
the strong clustering of haloes around the density peaks, the 
majority of these low-mass haloes are suppressed through- 
out most of the reionization process and hence their effect 
is mostly limited to the early stages of reionization and 
the progress of reionization is delayed. The time of over- 
lap and the large-scale geometry of reionization are both 
largely dominated by the high-mass haloes and thus are not 
strongly affected by the absence of low-mass haloes in the 
simulations. 

Could the inclusion of low-mass haloes and minihaloes 
significantly modify our results? A halo with a mass of 
10* M0 can reionize a mass between 2.5 x lO"^ -2.5 x IO^Mq 
(somewhat decreased once recombinations are accounted 
for) during its lifetime of ~20 Myr (before star formation 
shuts off again due to feedback, and recombinations set in, 
decreasing the reionized fraction further). Thus, in order to 
reionize a (proto-)halo with a mass around 1O^^M0, several 
10* Mq haloes ought to have formed at very similar times 
within the region considered. They should also have been 
spaced far enough apart so that the halo forming first does 
not suppress star formation in the other haloes. Even if this 
were the case, this group of ~1O*M0 haloes would most 
probably quickly merge, forming a halo with a mass above 
10^ Mq, which would be resolved in our simulations. More- 
over, regions which contain only 10* — 1O®M0 haloes, but not 
any of the more massive ones tend to be near the boundaries 
of the ionized regions, and are thus quickly overrun by the 
expanding I-fronts. Furthermore, both low-mass atomically- 
cooling haloes and minihaloes are strongly clustered around 
the density peaks, which are the first regions to be reionized, 
as we have shown. Therefore the vast majority of these small 
haloes would be affected by Jeans-mass filtering. As a conse- 
quence, haloes forming after their vicinity is ionized would 
not be able to accrete gas, while haloes that have formed 
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before their region became ionized would either not accrete 
more gas and thus have their star formation shut off (in the 
case of low-mass sources) , or would be photoevaporated (in 
the case of minihaloes). 

For all these reasons, we expect that neither haloes 
with masses between 10® — 10^ Mq, nor minihaloes would 
change our average reionization redshifts significantly and 
hence would not modify our results for the external/internal 
reionization fractions. Once future cosmological simulations 
become large enough to allow resolving low-mass haloes in 
a large volume, it would be interesting to verify if our as- 
sertions here are correct and evaluate the level of this effect 
precisely. 

Furthermore, the reionization parameters, photon pro- 
duction efficiencies and small-scale gas clumping (and to 
lesser level, also the background cosmology) are still quite 
uncertain. We have addressed this by comparing three 
reionization scenarios which roughly bracket the expected 
range of reionization histories, from the very early reion- 
ization given by f2000_WMAPl (high-efficiency sources, no 
sub-grid clumping, WMAPl cosmology), to the late over- 
lap, extended reionization scenarios, f250C_WMAPl and 
f250C_WMAP3 (low-efficiency sources, sub-grid clumping, 
WMAPl or WMAP3 cosmology). Most of the reasonable 
reionization scenarios which satisfy the current observa- 
tional constraints should fall between these cases. Further 
improvement would require better, more detailed observa- 
tional data which would constrain the reionization scenarios 
(and thus the reionization parameters) better. 



6 IMPLICATIONS OF EXTERNAL VS. 

INTERNAL REIONIZATION REGARDING 
GALAXY PROPERTIES 

Reionization can affect the properties of galaxies in many 
ways, and the relative timing of reionization across the mass 
scales will leave distinct observable signatures. 

• The mode of galaxy formation may be fundamentally 
different between haloes less massive and more massive than 
10^^ Mq. Late forming haloes that are externally ionized will 
have hotter, more extended and smoother gas distributions. 
The gas accretion into the galaxy will proceed differently, 
perhaps providing an explanation for the red-blue galaxy 
sequence, bulgeless disks, lower-baryon fractions in less mas- 
sive galaxies, etc. 

• Early star formation will be affected, since stars form- 
ing in ionized gas may have different p roperties than stars 
formin g in neutral gas. According to Ijohnson fc Bromml 
l|2006t ). Pop 11.5 could have formed from primordial gas that 
is strongly ionized since the free-electron-catalyzed forma- 
tion of molecules enables e nhanced gas cooling. Similarly, 
iNagakura fc Omukail l|2005l ) predict the formation of low 
mass stars in ionized gas with very low metallicity. External 
reionization provides an alternative pathway to Pop 11.5 star 
formation. 

• The abundance of Galactic satellites has been recon- 
ciled with observational data using the idea that reioniza- 
tion suppresses the formation of late forming haloes - only 
the rare peaks that collapse before reionization can cool gas 
to form stars. If the gas in a late-forming galaxy is externally 



reionized it may suppress the formation of some or even all 
of its satellite population. 

• The radial distribution of old stars in the Galac- 
tic Halo might refiect the rarity of the density peaks 
in which they formed at the epoch of reioniza tion 
(|Diemand. Madau. fc Moord |2005| . iMoore et all 120061 ). If 
stars form only in the rarest peaks, then their final distribu- 
tion can be highly concentrated with a radial profile as steep 
as r~*. Thus, externally reionized galaxies may be expected 
to have steeper diffuse stellar halo profiles. 

• Similar to the old stellar halo, globular cluster forma- 
tion may be truncated at reionization, altering the relative 
abundances and radial distributions of the old blue pop- 
ulation in g a laxies of different masses and environments. 
iRhode et at] (120051 ) found that the number of old blue GCs 
per galaxy stellar mass increases with increasing mass for 
cluster ellipticals, while it seems to be constant for spiral 
galaxies in the field. They suggest that the reason for this 
is that more massive elliptical galaxies assemble earlier and 
have more time to form blue GC before reionization, as- 
suming that reionization happens at a fixed point in time. 
Although we find that reionization will occur considerably 
earlier in massive elliptical galaxies than for low-mass field 
galaxies, their redshift of formation is shifted even more 
strongly, as they are usually internally reionized. Thus it 
is indeed the case that typically more time elapses between 
the formation and reionization of massive cluster ellipticals 
than for field galaxies. 

We conclude that it may already be possible to discrimi- 
nate between galaxies that have been reionized internally or 
externally using one or more of the above characteristics. 
This will be an interesting new and complementary way to 
probe the history of reionization in the Universe and its ef- 
fect on galaxy formation. 



7 SUMMARY AND DISCUSSION 

The reionization history of a given galaxy is strongly de- 
pendent on its mass due to the correlation between galaxy 
mass and formation redshift. Higher mass haloes on aver- 
age form earlier and are also reionized earlier. The lower the 
mass of a halo, the more important does external reioniza- 
tion (i.e. reionization emanating from sources outside the 
halo) become. The effect of the large-scale environment on 
the reionization history, on the other hand, is negligible for 
field haloes. By using a simple polynomial fit, we can pa- 
rameterize the average reionization redshift of a galaxy as a 
function of halo mass for three different combinations of cos- 
mology and reionizing photon production efficiency. Typical 
L* galaxies and LG like systems are reaching reionization 
considerably before the final overlap, at a redshift where 
roughly half of the Universe is reionized. Very massive, el- 
liptical galaxies in clusters are reionized even earlier and are 
always internally-reionized. Our prescription for the reion- 
ization redshift as a function of galaxy mass might be of use 
for semi-analytical models of galaxy formation. Normally, 
these models assume a fix ed, fairly low redshi ft of reioniza- 
tion of around z ~7 (e.g. ICroton et al.|[200^ . from which 
the baryon fraction in a low mass haloes can be calculated, 
given redshift and halo mass l|Gnedinll2000l . iKravtsov et al.l 
|2004 ). Using our fitting formulae, one could utilize a more 
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realistic reionization redshifts in semi-analytical models by 
introducing a dependence on the mass of the halo in which a 
given subhalo ends up. This may help alleviate the problem 
of the overproduction of sa tellite galaxies in massi ve clusters 
in semi-analytical models (IWeinmann et aLll2006l) . 

We suggest that whether a galaxy experiences internal 
or external reionization might affect its properties consider- 
ably. Externally reionized haloes could experience different 
gas accretion; could have a different first generation of stars 
(i.e. Pop II. 5 instead of Pop III) ; may host a lower number of 
satellite galaxies and blue globular clusters; and could have 
steeper diffuse stellar halo profiles. 
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